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Abstract. This paper is a continuation of our earlier paper in which we have 
derived the solution of an unified fractional reaction-diffusion equation asso- 
ciated with the Caputo derivative as the time-derivative and the Riesz-Feller 
fractional derivative as the space-derivative. In this paper, we consider an uni- 
fied reaction-diffusion equation with Riemann-Liouville fractional derivative as 
the time-derivative and Riesz-Feller derivative as the spacederivative. The so- 
lution is derived by the application of the Laplace and Fourier transforms in a 
compact and closed form in terms of the H-function. The results derived are 
of general character and include the results investigated earlier by Kilbas et al. 
(2006a), Saxena et al. (2006c), and Mathai et al. (2010). The main result is 
given in the form of a theorem. A number of interesting special cases of the 
theorem are also given as corollaries. 

1 Introduction 

In recent years, fractional reaction-diffusion models are studied due to their use- 
fulness and importance in many areas of mathematics, statistics, physics, and 
chemistry (Mainardi, 2010; Mathai and Haubold, 2008; Haubold and Mathai, 
2010; Mathai et al. 2010). Such models, formulated in standard or fractional 
calculus, greatly contribute to the understanding of the behavior of many-body 
systems far from equilibrium and emerging spatio-temporal pattern formation. 
Recently, coupled fractional reactiondiffusion equations are solved by Gafiychuk 
et al. (2006). Turing pattern formation through linear stability analysis and 
numerical simulation are discussed by Langlands et al. (2007). Nonlinear oscil- 
lations and stability domains in fractional diffusion systems with two types of 
variables activator and inhibitor are demonstrated by Gafiychuk et al. (2007). 
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General models for reaction-diffusion systems are investigated by Henry and 
Wearne(2000, 2002) and Henry ct al. (2005). 

This paper deals with the investigation of the solution of an unified model 
of reactiondiffusion system associated with the Riemann-Liouville fractional 
derivative as the time-derivative and the Riesz-Feller derivative as the space- 
derivative. This new model provides the extension of the models discussed ear- 
lier by Kilbas et al. (2006a), Saxena et al. (2006c), and Mathai et al. (2010). As 
special cases of this general model, we discuss neutral fractional diffusion model, 
time-fractional diffusion model and space-time fractional diffusion model. The 
present study is in continuation of our earlier works, Haubold and Mathai (f 995, 
2000, 2008, 20f0), Haubold et al. (2007, 20f0) and Saxena ct al. (2006a, 2006b, 
2006c). 

2 Results Required in the Sequel 

The RiemannLiouville fractional integral of order v is defined by Miller and 
Ross (f993, p. 45) and Kilbas et al. (2006) 

Dt v N(x, i) = J- /" (t - u)"- 1 ^*, u)du, (f ) 

1 \y) Jo 

where Re{v) > 0. 

The Riemann-Liouville fractional derivative of order a > is defined as 
(Samko et al., 1990, p.37; see also Kilbas et al., 2006) 

D?f{x,t) = — — f K x ^)dT = ^ 

tJK ' r(n- a) dt n J {t-T) a+1 - n w 

where [a] means the integral part of the number a. 

The Laplace transform of the Riemann-Liouville fractional derivative is given 
by Oldham and Spanier (1974, eq.(3.1.3); see also Kilbas et al., 2006) 

n 

L { D?N(x, t); s} = s a N(x, s) - ^ s^ 1 N(x, t)\t = 0(n - K a < n). 

r=l 

(3) 

The Riesz-Feller space-fractional derivative of order a and skewness is defined 
in terms of its Fourier transform as (Feller, 1952, 1971): 

F{ x D%f(x);k} = -$i(k)f*(k), (4) 

where f*(k) denotes the Fourier transform of /(t), defined by 

/oo 
exp(-ikt)f(t)dt (5) 
-oo 
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and 

^ e a (k) = \k\ a exp[i{signk)—],d<a<2, \6\ < min {a,2 - a} . (6) 

When 9 = 0, we have a symmetric operator with respect to x that can be 
interpreted as 

— (Hp) ■ m 

This can be formally deduced by writing — {k)a = — (fc 2 ) Q / 2 . Eq. (4) then 
reduces to 

F{ x Dgf(x);k} = -\k\ a , (8) 
which is the Fourier transform of the Weyl fractional operator, defined by 

p* m= 1 dU [* f^ du ( 9 ) 

This shows that Ricsz-Fcller operator may be regarded as a generalization of 
Weyl operator. An alternative notation for the symmetric fractional derivative 
$Dq has been given by Saichev and Zaslavsky (1997) in the form 

° D ° = < 10 » 

In its regularized form, which holds for < a < 2, the RieszFcllcr derivative 
admits the explicit representation 

, DSm . Hi±^,„ («) ^ /(x + a-2/( : ) + / ( x-o ^ 

For a = 1, the Riesz derivative is related to the Hilbert transform as pointed 
out by Feller (1952). We have 

7T dx J_ OQ x-£ 

For < a < 2 and |r/| < min {a, 2 — a}, the Riesz-Feller derivative can be 
shown to possess the following integral representation in x domain: 

xD?f(x) 

= {««[(« + *W2] ^ f{X + ® + ~ f{x) <% + sin[{a 9)n/2] 

f(x - - /(*) 



/•OO 



1+a 



de (13) 



3 



3 Unified Fractional Reaction-Diffusion Equa- 
tion 



In this section, we will investigate the solution of the reaction-diffusion system 

(14) under the initial conditions (15). 

Theorem. Consider the unified fractional reaction-diffusion model associated 
with Riemann-Liouville fractional derivative qD" defined by (2) and the Riesz- 
Feller space fractional derivative x Dg of order a and asymmetry 9 defined by 
(4) 

D?N(x, t) = ri x D$N(x, t) + $(x, t), (14) 
where r),t > 0,x € R;a,0, (3 are real parameters with the constraints 

< a < 2, \6\ < min(a, 2 - a), 1 < /3 < 2, 

and the initial conditions 

oD^Nix, 0) = f(x), D?~ 2 N(x, 0) = g{x) for x e R, lim N(x, t)=0,t>0. 

(15) 

Here Df _1 iV(a;, 0) means the Riemann-Liouville fractional partial derivative 
of N(x,t) with respect to t of order /? — 1 evaluated at t = 0. Similarly 
oD^ 2 N(x, 0) is the Riemann-Liouville fractional partial derivative of N(x, t) 
with respect to t of order j3 — 2 evaluated at t — 0. In (13), 77 is a diffusion 
constant and $(a:,t) is a nonlinear function belonging to the area of reaction- 
diffusion. Then for the solution of (14), subject to the above constraints, there 
holds the formula 

fP-i r°° 

N(x,t) = — / r(k)E M (- V t^ e a (k))exp(-ikx)dk (16) 
+ — J tg%k)E PiP - X (-rfi p il>i(k))exp(ikx)dk 

S^dt / <p*(k,t - C)E ^{-r l ^i } 6 a {k))exp{-ikx)dk, 

J — OO 

where E a ^(z) is the generalized Mittag-LefHer, defined by the series 

00 n 

E a ,p(z) = p7 " g) (a, (3) e C, Re(a) > 0, Re(J3) > 0. 

Proof. If we apply the Laplace transform with respect to the time variable t, 
Fourier transform with respect to space variable x and use the initial conditions 

(15) and the formula (17), then the given equation transforms into the form 

s^N? (k,s)- f* (k) - sg* (k) = -tjV£ (k)N* (k, s) + ^* (k,s), 
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where according to the conventions followed, the symbol ~ will stand for the 
Laplace transform with respect to time variable t and * represents the Fourier 
transform with respect to space variable x. Solving for N~(k,s) 7 it yields 

N *(k \= f * {k) I s 9*( k >^ 2 | ^( fc ) (m 
{,) sf> + V r a (k) sf> + 7?<(fc) + S P +rpl%(ky [ ' 

On taking the inverse Laplace transform of (16) by using the formula 

L ~" ( "4-77 1 = r- p Ea,a=fl+l{-<*r), (18) 



_ a + s' 

where Re(s) > 0,Re(a) > 0,Re(a — (3) > —1; it is seen that 

N *(k,t) = t^Vfi)^^-/ + i^V(fc)^,M(-//«W) 

Jo 

The required solution (16) is now obtained by taking the inverse Fourier trans- 
form of (19). This completes the proof of the theorem. 



4 Special Cases 

If we set 6 = by virtue of the results (6), the theorem reduces to the following 
Corollary 1. Consider the unified fractional reaction-diffusion model asso- 
ciated with Riemann-Liouville fractional derivative o-Df defined by (2) and the 
Riesz-Feller space fractional derivative x DqoI order a defined by (11) 

D?N(x, t) = Vx D°N(x, t) + <p(x, t), (20) 

where 77, t > 0, x E R; a, (3 arc real parameters with the constraints 

0<a<2,l</3<2, 



and the initial conditions 

o£>f ~ X N{x, 0) = f(x), D f t~ 2 N(x, 0) = g(x) iorx e R, lim N(x, t)=0,t>0. 

(21) 

Here Q D% ~ X N{x, 0) means the Riemann-Liouville fractional partial derivative 
of N(x,t) with respect to t of order (3 — 1 evaluated at t = 0. Similarly 
oDf~ 2 N(x, 0) is the Riemann-Liouville fractional partial derivative of N(x,t) 
with respect to t of order (3 — 2 evaluated at t = 0. The r\ is a diffusion constant 
and <fi(x,t) is a nonlinear function belonging to the area of reaction-diffusion. 
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Then for the solution of (14), subject to the above constraints, there holds the 
formula 



N(x,t) 



+ 



2tt 

t f)-2 



2ir 



1 roo 

- I f*{k)E (j j j {-r 1 t fi \k\ a )exp{-ikx)dk 

J — OO 

o 

g*(k)E ^_ 1 (-r]t l3 \k\ a )exp(-ikx)dk 



+ 2^J J ^( k ' t -OE A-V^\k\ a )exp(-ikx)dJ^2) 

When g(x) = 0, then by the application of the convolution theorem of the 
Fourier transform to the solution (16) of the theorem, it readily gives 

Corollary 2. The solution of fractional reactiondiffusion equation 

D^N(x, t)-r) x D%N(x, t) = ip(x, t), x G R, t > 0, rj > 0, (23) 

with initial conditions 

oDf^NfaO) = f(x), Dt 2 N(x,0) - 0, 

for x G R, - < a < 1, 1< /3 < 2, lim iV(a;,i)=0, (24) 

where r\ is a diffusion constant and (fi(x,t) is a nonlinear function belonging to 
the area of reaction-diffusion; r),t > 0,x G R; a, 0,(3 are real parameters with 
the constraints 



is given by 

N(x,t) 



< a < 2, \0\ < min{a, 2 - a), 1< /3 < 2 

= / Gi (x — t, t)f(r)dT 
Jo 

+ f dtit-Zf- 1 ^ G 2 (x-T,t-0<p(T,0dT, 

Jo Jo 



(25) 



where 



P 

G!(x,t) 



a-6 

tP- 1 
-rrH. 



exp{~ikx)E fi ^{-rit f} il) f> a {k))dk 



a\x\ 



2,1 
3,3 



(l,l/a),(/3,/3/a),(l,p) 



, (a > 0) (26) 
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and 



i r°° 

G 2 {x,t) = — / exp{-ikx)Ep, fj {~r ] t^ e a (k))dk 

IK J-oo 



a\x 



1 w 2 ' 1 

— ^3,3 



v l/ at P/c 



(l,l/a),(/3,0/a),(l,p) 



(a > 0). (27) 



Here iJ 3 ' 3 (.) is the H-function , denned in Mathai and Saxena (1978, p. 2; see 
also Mathai et al. 2010). In deriving the above results, we have used the inverse 
Fourier transform formula (Haubold et al., 2007) 



F- 1 [E^{-^{k));xj\ = 4^H 2 3: l 



(l,l/a),( 7 ,/3/a),(l,p) 
(l,l/a),(l,l),(l,p) 



(28) 



where Re(a) > Re{/3) > 0, Re{i) > 0. 



It is interesting to observe that for 9 = Corollary 1 reduces to a result given 
by the authors (Saxena et al., 2006c). On the other hand if we set f(x) = (r(x), 
where a(x) is the Dirac-delta function, it yields 

Corollary 3. Consider the following reaction-diffusion model 

D?N(x,t)=rixD?N(x,t), (29) 

with the initial conditions n D^ 1 N(x, 0) = a(x),0 < (3 < 1, lim^-too N(x, t) = 
0, where rj is a diffusion constant; r],t > 0,x G R;a 7 8 7 (3 are real parameters 
with the constraints 



< a < 2, \6\ < min(a,2- a) 



and 6(x) is the Dirac-delta function. Then for the fundamental solution of (29) 
with initial conditions, there holds the formula 



2,1 



(^/3)l/« 



(l,l/a),(l,l),(l,p) 



, (a > 0) (30) 



where p = 3 y^-- For 9 = we obtain the result given by Kilbas et al. (2006a) 
obtained in a different form. In this case the authors have given the result in a 
closed form (Saxena et al., 2006c, p. 309). 

Remark. We note that the equation (29), when the Riemann-Liouville operator 

appearing on its left is replaced by a Caputo derivative (Caputo, 1969; Mainardi, 

2010), has been solved by Mainardi et al. (2001, 2005). 

The following special cases of (29) are worth mentioning: 

(i) For a = /3, the corresponding solution of (29), denoted by , we call 

as the neutral fractional diffusion, which can be expressed in terms of the H- 

function as given below and can be defined for x > 0: 
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Neutral fractional diffusion: 0<a~f3<2;9< min {a, 2 — a} , 

t a-l 



-H. 



2.1 

3,3 



^1/alpha 



(l,l/cO,(a,l),(l,p) 
(l,l/a),(l,l),(l,p) 



P = 



2a 



(31) 



Next we derive some stable densities in terms of the H-function as special 
cases of the solution of the equation (26). (ii) When (3 = 1,0 < a < 2;9 < 
min {a, 2 — ajthen (29) reduces to space-fractional diffusion equation, which 
we denote by is the fundamental solution of the following space-time fractional 
diffusion model: 



dN(x,t) 
dt 



r] x D%N(x,t),ri>0,x€R 



(32) 



with the initial conditions N(x, t = 0) = (t(x) 7 \\ir x ^± 00 N(x, t) = 0, where rj is 
a diffusion constant and 5(x) is the Dirac-delta function. Hence for the solution 
of (29) there holds the formula 



1 _ TJ-1,1 

a{vt) l/ a ^2,2 



{fft) 



1/a 



(1,1),(P.P) 



< a < 1, \0\ < a, (33) 



where p = S 2a~- The density represented by the above expression is known as 
a-stable Levy density. Another form of this density is given by 



L 6 a (x) = 



a{r]t) l/a H 2.2 



(0,l),(l-p,p) 



, 1< a2, \6\ < 2 - a. (34) 



Note: A comprehensive account of stable densities with applications is available 
from the monograph of Uchaikin and Zolotarev (1999). 

(iii) Next, if we take a = 2, < (3 < 2; 9 = 0, then we obtain the time- 
fractional diffusion, which is governed by the following time-fractional diffusion 
model: 

d fi N( x ,t) d 2 



ri—N{x,t),f] > 0,x e R, < /3 < 2, 



(35) 



with the initial conditions o-Df ~ 1 N(x, 0) = a(x),o £>f ~ 2 N(x, 0) = 0, for x G 
r, lim^^ioo N(x, t) — where eta is a diffusion constant and a(x) is the Dirac- 
delta function, whose fundamental solution is given by the equation 



N(x,t) 



2\x\ 



tt\,0 



(W(l,l) 



(36) 



Further, if we set a — 2, (3 = 1 and 9 — > 0, then for the fundamental solution of 
the standard diffusion equation 



^N { x, t )=n^N(x,t), 



(37) 
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with initial condition 



N(x, t = 0)= a(x), lim N(x, t) = 



(38) 



there holds the formula 



N{x 1 t) = —H{f 



^1/2^/2 (1,1) 



N (1,1/2) 



(47T?7£) 1/2 ea;p[ 




(39) 



which is the classical Gaussian density. 

In conclusion, it is seen that the solution given by (27) does not admit a 
probabilistic interpretation in contrast with fractional reaction-diffusion based 
on Caputo derivative derived by the authors (Haubold et al., 2007). However, 
when (3 — > 1 , then it has a probabilistic interpretation, as can be seen in special 
cases of corollary 3. 
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